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Abstract. In this paper, we consider the wave equation with both a viscous Kelvin- Voigt and 
frictional damping as a model of viscoelasticity in which we incorporate an internal control with 
a moving support. We prove the null controllability when the control region, driven by the flow 
of an ODE, covers all the domain. The proof is based upon the interpretation of the system 
as, roughly, the coupling of a heat equation with an ordinary differential equation (ODE). The 
presence of the ODE for which there is no propagation along the space variable makes the 
controllability of the system impossible when the control is confined into a subset in space that 
does not move. The null controllability of the system with a moving control is established in 
using the observability of the adjoint system and some Carleman estimates for a coupled system 
of a parabolic equation and an ODE with the same singular weight, adapted to the geometry 
of the moving support of the control. This extends to the multi-dimensional case the results by 
P. Martin et al. on the one-dimensional case, employing 1 — d Fourier analysis techniques. 



1. Introduction 

We are concerned with the controllability of the following model of viscoelasticity consisting 
^ of a wave equation with both viscous Kelvin- Voigt and frictional damping: 

ytt ~ Ay - Ay t + b{x)y t = l w(t) /i, x€£l,t€(0,T), (1.1) 
l-Ej y = 0, xedQ, t€(0,T), (1.2) 

i— h y(x,Q) = y (x), y t (x,0) = yi(x), x e fl (1.3) 

^ Here Q, is a smooth, bounded open set in R N , b G L°°(Q) is a given function determining 
the frictional damping and h = h(x, t) denotes the control. To simplify the presentation and 
notation, and without loss of generality, the viscous constant has been taken to be the unit one 

^] v = 1. The same system could be considered with an arbitrary viscosity constant u > leading 

^ to the more general system 

cn yu - Ay - uAy t + b(x)y t = l^h, (1.4) 

^ but the analysis would be the same. 

The control h acting on the right hand side term as an external force is, for all < t < T, 
localized in a subset of £1. This fact is modeled by the multiplicative factor l w (t) which stands 
for the characteristic function of the set u(t) that, for any < t < T, constitutes the support of 
the control, localized in a moving subset u(t) of U. 
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Typically we shall consider control sets uj(t) determined by the evolution of a given reference 
subset u of Q through a smooth flow X(x, t, 0). 

We consider the problem of null controllability. In other words, given a final time T and 
initial data for the system (yo,Ui) in a suitable functional setting, we analyze the existence of a 
control h = h(x,t) such that the corresponding solution satisfies the rest condition at the final 
time t = T: 

y(x,T) = y t (x,T) = 0, in O. 

One of the distinguished features of the system under consideration is that, for this null 
controllability condition to be fulfilled, the control needs to move in time. Indeed, if ui(t) = uj 
for all < t < T, i.e. if the support of the control does not move in time as it is often 
considered, the system under consideration is not controllable. This can be easily seen at the 
level of the dual observability problem. In fact, the structure of the underlying PDE operator 
and, in particular, the existence of time-like characteristic hyperplanes, makes impossible the 
propagation of information in the space-like directions, thus making the observability inequality 
also impossible. This was already observed in the work by P. Martin et al. in [ ] in the 1 — d 
setting. There, for the 1 — d model, it was shown that this obstruction could be removed by 
making the control move so that its support covers the whole domain where the equation evolves. 

More precisely, in [25], the 1 — d version of the problem above was considered in the torus, 
with periodic boundary conditions, 6 = and oj{t) = {x — t; x £ u}, i.e. 

Vtt ~ Vxx ~ Vxxt = l^(t)h{x, t), x e T. (1.5) 

Recall that this system with boundary control, i.e. h = and the boundary conditions 

y(0,t) = 0, y(l,t)=g(t), 

g = g(t) being the boundary control, fails to be spectrally controllable, because of the existence 
of a limit point in the spectrum of the adjoint system [28]. In the moving frame x' = x + t, (1.5) 
may be written as 

z tt - 2z xt - z txx + z xxx = a(x)h(x + t, t) (1.6) 
where z(x, t) = y(x + t,t). In [25] the spectrum of the adjoint system to (1.6) was shown to be 
split into a hyperbolic part and a parabolic one. As a consequence, equation (1.6) was proved 
to be null controllable in large time. A similar result was proved in [31] for the Benjamin-Bona- 
Mahony equation 

yt - ytxx + y x + yy x = a(x - ct)h(x, t), x g T. 
Once again this system turns out to be globally controllable and exponentially stabilizable in 
ff 1 (T) for any c / 0. But, as noticed in [26], the linearized equation fails to be spectrally 
controllable with a control supported in a fixed domain. 

As mentioned above, in both cases, the lack of controllability of these systems with immobile 
controls is due to the fact that the underlying PDE operators exhibit the presence of time-like 
characteristic lines thus making propagation in the space-like directions impossible. By the 
contrary, when analyzing the problem in a moving frame, the characteristic lines are oblique 
ones in (x,t), thus facilitating propagation properties. 

The main goal of this paper is to extend the 1 — d analysis in [25] to the multi-dimensional 
case. This can not be done with the techniques in [25] based on Fourier analysis. Our approach 
is rather inspired on the fact that system (1.1)-(1.2) can be rewritten as a system coupling a 
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parabolic equation with an ordinary differential equation (ODE). The presence of this ODE, 
in the case of a fixed support of the control, independent of t, is responsible for the lack of 
controllability of the system, due to the absence of propagation in the space-like direction. 
Letting the control move introduces an effect similar to adding a transport term in the ODE 
but keeping the control immobile, thus changing the structure of the system into a parabolic- 
transport coupled one. This new system turns out to be controllable under the condition that all 
characteristics of the transport equation enter within the control set in the given control time, 
a condition that is reminiscent of the so-called Geometric Control Condition in the context of 
the wave equation (see [2]). 

The approach in [25] would suggest to do the following splitting of (1.1): 

vt-Av = i u{t) h + (l-b)(v-y) (1.7) 

y t + y = v. (1.8) 

However, the splitting can be performed in an alternative manner as follows: 

y t - Ay + (6 - l)y = z (1.9) 

zt + z = l u(t) h + (b-l)y. (1.10) 

It can easily be seen that y solves equation (1.1) if, and only if, it is the first component of the 
solution of system (1.9)-(1.10). 

Our analysis of the Carleman inequalities for these systems is analog to that in [ ] for a 
system of thermoelasticity coupling the heat and the wave equation. The key in [1] and in our 
own analysis is to use the same weight function both for the Carleman inequality of the heat 
and the hyperbolic model. In [1], since dealing with the wave equation, rather strong geometric 
conditions were needed on the subset where the control or the observation mechanism acts. In 
our case, since we are considering the simpler transport equation, the geometric assumptions will 
be milder, consisting mainly on a characteristic condition ensuring that all characteristic lines 
of the transport equation intersect the control/observation set. This suffices for the Carleman 
inequality to hold for the transport equation and is also sufficient for the heat equation that it 
is well-known to be controllable/observable from any open non-empty subset of the space-time 
cylinder where the equation is formulated. 

It is important to mention that, as far as we know, all the Carleman inequalities for the heat 
equation available in the literature are done for the case where the control region is fixed. In the 
case we are dealing, the control region is moving in time. Therefore, a new Carleman inequality 
must be proved in this framework. The proof of a Carleman inequality for the heat equation 
when the control region is moving is one of the novelties of this paper. 

In order to state the main result of this paper we first describe precisely the class of moving 
trajectories for the control for which our null controllability result will hold. 

Admissible trajectories: In practice the trajectory of the control can be taken to be 
determined by the flow X(x,t,t ) generated by some vector field / G C([0,T]; W 2 '°°(R N ;R N )), 
i.e. X solves 

dX 

— (x, t, t ) = f(X(x, t, to),t), 
X(x,t ,t ) = x. 
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For instance, any translation of the form: 

X(x,t,t ) = x + j(t) -j(t ), (1.12) 

where 7 G C 1 ([0, T]; M^), is admissible. (Pick = 7(4)). 

We assume that there exist a bounded, smooth, open set cj C K^, a curve T G C°°([0, T]; R N ), 
and two times ti, t<i with < ti < ti < T such that: 

r(t)e%i,o)nfi, vte[o,r]; (1.13) 

O C U mo>T] X(co , t, 0) = {X{x, t, 0); x G u , t G [0, T]}; (1.14) 



17 \ X(wo, i, 0) is nonempty and connected for t G [0, ti] U [t2, T\; (1-15) 

£1 \ X(ujo,t,0) has two (nonempty) connected components for t G (ti, ^2); (1-16) 

V7GC([0,r];J2), 3iG [0,T], 7 (t) G X(a; ,i,0). (1.17) 
The main result in this paper is as follows. 

Theorem 1.1. Le£ T > 0, X(a;,i,io) ^0 be as in (1.13)-(1.17), and let oj be any open set 
in Q, such that loo C u. Then for all (yo,yi) G L 2 (Q,) 2 with y\ — Ayo G L 2 (Q), there exists a 
function h G L 2 (0,T; L 2 (Q)) for which the solution of 

y tt - Ay- Ay t + b(x)y t = l x{ujm ( y x)h, (x, t) G Q X (0, T), (1.18) 

y(x, t) = 0, (x, t)£dttx (0, T), (1.19) 

y(.,0)=y , y t (.,0) = yi , (1.20) 

fulfills y(.,T) = y t (.,T) = 0. 

A few remarks are in order in what concerns the functional setting of this model: 

• Viewing the system of viscoelasticity under consideration as a damped wave equation, 
a natural functional setting would be the following: For data in Aq(O) x L 2 {yi) and, 
say, right hand side term of (1.1) in L 2 (0,T; H^ 1 there exists an unique solution 
y G C([0,T];H^(Q) n C x ([0, T]; L 2 (Q)). Furthermore y t G L 2 (0, T; ^(0)). The latter 
is an added integrability/regularity property of the solution that is due to the strong 
damping effect of the system. This can be seen naturally by considering the energy of 
the system 

E(t)= 1 -Jj\y t \ 2 + \Vy\ 2 ]dx, 



that fulfills the energy dissipation law 
dt 



E(t) = - f [\Vy t \ 2 + b(x)\y t \ 2 ]dx + / hy t dx. 
Jn Ju){t) 



We can also solve (1.9)-(1.10) so that y, solution of the heat equation, lies in the space 
y G C([0,T];H^(n))nL 2 (0,T;H 2 (n)) and z, solution of the ODE, in C([0, T\; L 2 (Q)). 
This can be done provided (yo,yi — Ayo) G Hq(CI) x L 2 (Q). The functional setting is 
not exactly the same but this is due to the fact that, in some sense, in one case we see 
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the system as a perturbation of the wave equation, while, in the other one, as a variant 
of the heat equation. 

• From a control theoretical point of view it is much more efficient to analyze the system 
in the second setting, as a perturbation of the heat equation, through the coupling with 
the ODE or, after changing variables, with a transport equation. If we view the system 
of viscoelasticity as a perturbation of the wave equation, standard hyperbolic techniques 
such as multiplier, Carleman inequalities or microlocal tools do not apply since, actually, 
the viscoelastic term determines the principal part of the underlying PDE operator and 
cannot be viewed as a perturbation of the wave dynamics. 

The analysis is particularly simple in the special case where 6 = 1. Indeed, in that case, the 
system (1.7)-(1.8) (with a second control incorporated in the ODE) takes the following cascade 
form 

v t -Av = l u(t) h, (1.21) 

yt + y = i u (t)k + v, (1.22) 

where the parabolic equation (1.21) is uncoupled. This system will be investigated in a separate 
section (Section 2) since some of the basic ideas allowing to handle the general case emerge 
already in its analysis. Note that, in this particular case, roughly, one can first control the heat 
equation by a suitable control h and then, once this is done, viewing v as a given source term, 
control the transport equation by a convenient k. This case is also important because the only 
assumption needed to prove Theorem 1.1 in this case is (1-14) (i.e. we don't assume that (1.13) 
and (1.15)-(1.17) are satisfied). 

In this particular case b = 1 a similar argument can be used with the second decomposition. 

The paper is organized as follows. Section 2 is devoted to address the particular case 6=1. 
In Section 3 we give some examples showing the importance of the taken assumptions on the 
trajectories. In Section 4, we go back to the general system (1.9)-(1.10). We prove that this 
system is null controllable in L 2 (Cl) 2 (see Theorem 4.1) by deriving the observability of the 
adjoint system from two Carleman estimates with the same singular weight, adapted to the flow 
determining the moving control. The details of the construction of the weight function based on 
(1.13)-(1.17) are given in Lemma 4.3. Theorem 1.1 is then a direct consequence of Theorem 4.1. 
We finish this paper with two further sections devoted to comment some closely related issues 
and open problems 

2. Analysis of the decoupled cascade system 

In this section, we give a proof of Theorem 1.1 in the special situation when 6 = 1, and 
w(t) = X(uo,t,0), where X is given by (1.11) for some / G C(R+; W 2 >°°(R N ; R N )). 

As we said before, we will prove Theorem 1.1 in the case 6 = 1 by proving a null controllability 
result for the decomposition (1.21)-(1.22). The idea of the proof is as follows. We take some 
appropriate < e < T and then drive the solution of the heat equation (1.21) to zero in time 
e by means of a control h. Next, we let equation (1.22) evolves freely in [0, e], i.e., k = 0, and 
then we control this equation by means of a smooth control k in the time interval [e, T}. This 
gives the null controllability of the system (1.21)-(1.22) in the whole time interval [0, T]. 
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Proof of Theorem 1.1 in the case 6=1. 

Suppose (1.14) is satisfied and let 

T = inf{T > 0; ft C U < t <T^(w , *, 0)}. (2.1) 

Pick any T > To, and pick some e G (0, T — To) and some nonempty open set U-i C ujq such 
that 

ft C U e < t < T X(w ,i,0), (2.2) 

w_i C X(w ,t,0) Vt€(0,e). (2.3) 

Let T' = T — e, and pick any (vo,yo) G T 2 (ft) 2 . Then, it is well known (see [16]) that there 

exists some control input h G T 2 (0, e; L 2 (ft)) such that the solution v = v(x,t) of 

v t -Av = l u _ 1 h 1 x G ft, t G (0,e), (2.4) 

u(x,0) = «o(a0, ieO, (2.5) 

satisfies 

f (x, e) = 0, x G ft. 

Set 

h(x,t) = l LU _ 1 (x)h(x,t), t G (0, e), 

/i(x,t) = 0, x G ft, i G (e,T), 
k(x,t) = 0, x G ft, t G (0,e) 

Then the solution v of 

«t - At; = l X (u ,t,o)( x )h x G ft, t G (0,T), (2.6) 

u(o?,0) = «o(»)) x G ft, (2.7) 

satisfies v(x,t) = for i G [e, T]. We claim that the system 

ife + y = l x(oJom (x)k(x, t), x G ft, t G (e, T), (2.8) 

y(x,e) = y (x), (2.9) 

is exactly controllable in L 2 (ft) on the time interval (s,T). By duality, this is equivalent to 
proving that the corresponding observability inequality 

\q (x)\ 2 dx<C [ [ l x{uom (x) \q(x,t)\ 2 dxdt (2.10) 
Je Jn 

is fulfilled with a uniform constant C > for all solution q of the adjoint system 

-q t + q = 0, xGft, tG(e,T), (2.11) 

q(x,T)=q (x). (2.12) 

Since the solution of (2.11)-(2.12) is given by q(x,t) = e t ~ T qo(x), we have that 

f T f lx(u,o,t,o)(aO \q(x,t)\ 2 dxdt > e 2 ( £ - T ) / \q (x)\ 2 ( f lx (w0)t) o)(x) dt)dx. (2.13) 
Je Jfl Jfl Je 
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From (2.1) and the smoothness of X, we see that for all x G f2, there is some to £ ( e )^); 
and some 5 > such that for any y G B(x,5) and any t G (e, T) n (to — <5, to + <$) we have 
y G X(cJO) t, 0). From the compactness of f2, we see that there exists a number <5q > such that 

T _ 

1 X(w ,t,o)(a:)rf* > <5o, Vx G 0. 

Combined with (2.13), this yields (2.10). Thus, (2.8)-(2.9) is exactly controllable in L 2 (9) on 
(e,T) with some controls k G C([e,T]; L 2 (Q)). Let j/i(z) = e~ e yo(x) + e s_e f (x, s)e£s. Extend 
fc to (0,T) so that fc G L 2 (0, T; L 2 (0)) and the solution of 

Vt + y = lx(fj a ,t,o)(x)k, xen,te(e,T), (2.14) 
2/(x,e) = yi(x), (2.15) 

satisfies y(-,T) = 0. Thus the control (/i, k) steers the solution of (1.21)-(1.22) from (vo,yo) at 
t = to (0, 0) at t = T. Applying the operator dt — A in each side of (1.22) results in 

y tt -Ay- Ay t + y t = lx(u ,t,o)(x)h + {dt - A)[l x(uj(htfi) (x)k}, (2.16) 

This proves Theorem 1.1, except the fact that the control does not live in L 2 (0, T; L 2 (Q,)). 
Assume now that (v ,yo) G L 2 (Q) x [H 2 (Q) n H^(Q,)]. To get a control k G L 2 (0,T; H 2 (n)), it 
is sufficient to replace l w (t) by a(X(x,0,t)) in (1.22), where a is a function satisfying 

a(x) = 1 Vx G u>o- 
The proof is completed by showing the observability inequality 

WqoWh < C £ \\(t-e)a(X(.,0,t))q(.,t)\\ 2 x ,dt 

for the solution q of system (2.11)-(2.12), where X = H 2 {Q) n Hq(Q) and A" stands for its 
dual space. This can be done as in [ : , Proposition 2.1]. Next, using the HUM operator, we 
notice that k G C x {[e,T\;X) with k(.,e) = 0, since q G C x {[e,T\;X') for any q G X' . Thus, 
with this small change, the right hand side term in (2.16) can be written lx(w,t,o)' u ( x ) t), where 
u G L 2 (0,T;L 2 (Q)). 

Remark 2.1. Observe that the situation when X{uo,t, 0) moves as in Figure 2, Figure 4 or in 
Figure 5 (see below) is admissible in the case when 6=1. 

3. Examples 

In this section, we provide some geometric examples to illustrate the assumptions (1.13)- 
(1.17). We use simple shapes (like rectangles) just for convenience. 

• Figure 1 shows how a control region should move in order to satisfy conditions (1.13)- 
(1.17). 

• Figure 2 depicts a situation for which Theorem 1.1 cannot be applied, except in the case 
when b = 1, as condition (1.16) fails. 



X(« ,i,0) 
r >. 
r 

• r(t) 
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X(w o ,t,0) 



nj(i) 



r(i) 



< i < d 



t 1 <t <t 2 



t 2 <t<T 



Figure 1. Example for which conditions (1.13)-(1.17) are satisfied. 



Figure 2. Example for which condition (1.16) fails. 





XK*,o) 

Figure 3. Example for which condition (1.15) fails. 



X(uo,T,0) 




• In Figure 3, we modify the example given in Figure 1 by shifting the time. Theorem 
1.1 cannot be applied as it is, since (1.15) fails. However, the conclusion of Theorem 1.1 
remains valid. Indeed, assume that Q\u)(t) has two connected components (resp. one) 
for t £ [0, ti) U (t2,T] (resp. for t £ [ti,T2]). Assume that the "jump" of u)(t) occurs at 
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Figure 4. Example showing that X(uq, t, 0) n Q ^ Vi € [0, T] does not imply (1.13). 

t = T3, with n < T3 < T2- Let 

01 :=U <K T3 w(t), (3.1) 

2 := U r3 < t < T w(t) (3.2) 

and let 77 G C°°(ft; [0, 1]) be such that 

supp(r ? ) C Oi, (3.3) 
supp(l - t?) C 2 , (3.4) 
supp(Vr/) C ujq. (3.5) 

Then, applying the Carleman estimate in Lemma 4.6 to (pi,qi) = r)(X(x,0,t))(p,q) in 
0, n {r; > 0} on the time interval [0, T3], and to (^2,^2) = (l — f](X(x, 0, t))) (p, g) in 
Onjij < 1} on the time interval [T3, T], we can easily prove the observability inequality 
(4.11). 

• Figure 4 shows that the assumption (1.13), which is needed to construct the weight 
function tp in Lemma 4.3 cannot be replaced by the simpler condition 

i(w ,t,o)nfi^, Vie[o,T]. 

• Figure 5 shows that the assumption (1.17), which is also needed to construct the weight 
function iji in Lemma 4.3, does not result from the other assumptions (1.13)-(1.16). 

4. Null controllability of system (1.21)-(1.22). 

In this section we proof Theorem 1.1. Using decomposition (1.9)-(1.10), it is easy to see that 
the null controllability of (1.1)-(1.3) turns out to be equivalent to the null controllability of the 

system 
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X{u o ,t,0) X(u ,t,0) 




Figure 5. Example showing that (1.13)-(1.16) does not imply (1.17). 



y t - Ay + (b(x) - l)y = z, (x,t) € SI X (0,T), (4.1) 

zt + z = l x(ujm (x)h+(b(x)-l)y, Mefix(0,n (4.2) 

y(x,t) = 0, (x,t) € dSl X (0,T), (4.3) 

z(x,0) = afo(x), (4.4) 

y(x,0) = yoCa), (4.5) 
More precisely, Theorem 1.1 is a direct consequence of the following result. 

Theorem 4.1. LetT, X(x,t,to) and ojq be as in (1.13)-(1.17), and let u> be as in Theorem 1.1. 
Then for all (yo,Zo) 6 L 2 (il) 2 , i/iere exists a control function h € L 2 (0,T; L 2 (Sl)) for which the 
solution (y,z) of (4.1)-(4.5) satisfies y(.,T) = z(.,T) = 0. 



From now on we concentrate in the proof of Theorem 4.1. 
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It is well-known (see [12] ) that Theorem 4.1 is equivalent to prove an observability inequality 
for the adjoint system of (4.1)-(4.5), namely 

- Pt - Ap + (b(x) - l)p = (b(x)-l)q, (x,t) € O x (0,T), (4.6) 

-q t + q = p, (x,t)eOx(0,T), (4.7) 

p(x,t) = 0, (x,t) € dn x (0,T), (4.8) 

pfoT) = po(x), xeQ, (4.9) 

g(a;,r) = q (x), x £ O. (4.10) 

In fact, one can show that Theorem 4.1 is equivalent to the following: 

Proposition 4.2. Let T, X, ujq and uj be as in Theorem 1.1. Then there exists a constant 
C > such that for all (po,qo) £ L 2 (ft) 2 , the solution (p,q) of (4.6)-(4.10) satisfies 

[ [\p(x,0)\ 2 + \q(x,0)\ 2 ]dx < C [ [ \q(x,t)\ 2 dxdt. (4.11) 

Jn Jo Jx(u),t,o) 

Proof of Proposition 4-2. Inspired in part by [1] (which was concerned with a heat-wave system 1 ), 
we shall establish some Carleman estimates for the (backward) parabolic equation (4.6) and the 
ODE (4.7) with the same singular weight. 

For a better comprehension, the proof will be divided into two steps as follows: 

Step 1. We apply suitable Carleman estimates for the parabolic equation (4.6) and the ODE 
(4.7), with the same weights and with a moving control region. 

Step 2. We estimate a local integral of p in terms of a local integral of q and some small order 
terms. Finally, we combine all the estimates obtained in the first step and derive the desired 
Carleman inequality. 

The basic weight function we need in order to prove such inequalities is given by the following 
Lemma. 

Lemma 4.3. Let X, ujq and uj be as in Theorem 1.1, and let ui\ be a nonempty open set in M. N 
such that 

wq C u)\, U\ C co. (4-12) 



^See also [9] for some Carleman estimates for a coupled system of parabolic-hyperbolic equations. 
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t 

{(x,t); < t < T, x e X(u u t,Q)} 

T 



tpt<0 



A > 



Figure 6. Sign of the time derivative of ijj. 

Then there exist a number 5 £ (0, T/2) and a function ip £ C°°(f2 x [0, T]) such that 

Vi/>(x,t) / 0, t £ [0,T], x EQ\X(uji,t,0), 



ip(x,t) > 



ip t (x,t) > 0, 
ilh(x,t) < 0, 

^(M)<o, 

i°°(nx(o,T))) 



te [0,T], x £0\ A(wi,t,0), 
t £ [0,5], x £ H\ X(wi,t,0), 
te [T-<5,T], x £ n\X(wi,t,0), 

t £ [0,T], x £ <9ft, 
t £ [0,T], x £ 0. 



The proof of Lemma 4.3 will be given in Appendix A. 
Next, we pick a function g £ C°°(0,T) such that 



9(t) 



( l 

t 

strictly decreasing 
1 



and define the weights 



<p(x,t) 
8{x,t) 



{ giT-t) 



5 ( t )( e |A||^|U= 
g(t)e x ^ l \ 



for < t < 5/2, 
for < t < 5, 
for 5 < t < f , 
for £ < t < T 



(x,t) £ O x (0,T), 
(x,t) £ O x (0,T), 



where 



l°°(Qx(o,t)) an( i A > is a parameter. 



(4.13) 
(4.14) 
(4.15) 
(4.16) 

(4.17) 
(4.18) 
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Step 1. Carleman ESTIMATES with the same weight. 

In this step we apply a Carleman inequality for the heat-like equation (4.6) and a Carleman 
inequality for the ODE (4.7), both with the same weight. We combine such inequalities and 
obtain a global estimation of p and q in terms of local integrals of p and q. 

For the purpose of the proof, we assume that the following two lemmas are true (their proof 
are given, respectively, in Appendices B and C). 

Lemma 4.4. There exist some constants Xq > 0, so > and Cq > such that for all A > Ao, 
alls > s and all p G C([0, T]; L 2 (Q)) with p t + Ap G L 2 (0, T; L 2 (Q)), the following holds 

f T f [(sO)-\\Ap\ 2 + \ Pt \ 2 ) + A 2 (s^)|Vp| 2 + X 4 (s9) 3 \p\ 2 ]e- 2sip dxdt 
Jo Jn 

<C ([ T [\pt + Ap\ 2 e- 2s ^dxdt+[ T [ \ A {s6) 3 \p\ 2 e- 2s ^dxdt j . (4.19) 

\Jo Jn Jo Jx(uii,t,o) ) 



Lemma 4.5. There exist some numbers \\ > Ao, s\ > so and C\ > such that for all A > Ai, 
alls > si and all q G H 1 ^, T; L 2 (Q)), the following holds 

f T f (X 2 se)\q\ 2 e- 2sip dxdt<C 1 ( f T f \q t \ 2 e~ 2sLp dxdt+ F I \ 2 (s6) 2 \q\ 2 e - 2sip dxdt) . 

Jo Jn \Jo Jn Jo Jx(u,t,o) J 

(4.20) 

Applying the Carleman inequality given in Lemma 4.4 to the heat-like equation (4.6), we 
obtain 



[ T f [(sO)-H\^P\ 2 + \Pt\ 2 ) + X 2 (sO)\Vp\ 2 + \\s6f\p\ 2 ]e- 2 ^dxdt 
Jo Jn 



<C ( f f \{b{x-l){p-q)\ 2 e- 2sip dxdt+ T f \ A (s6) 3 \p\ 2 e- 2sip dxdt) . (4.21) 

Wo Jn Jo Jx(ui 1 ,t,o) I 



Next, we apply the Carleman inequality given by Lemma 4.5 to the ODE (4.7) and obtain 



(\ 2 s9)\q\ 2 e- 2s,f dxdt 



< Ci ( Ft \q-p\ 2 e- 2sip dxdt+ Ft \ 2 {se) 2 \q\ 2 e- 2s *dxdt j . (4.22) 

\Jo Jn Jo Jx(u,t,o) J 

Adding (4.21) and (4.22), it is not difficult to see that 
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Ff [(se)- 1 {\A P \ 2 + \p t \ 2 ) + \ 2 {s9)\V P \ 2 + \ 4 (sef\p\ 2 }e- 2s ^dxdt+ Ff {\ 2 s9)\q\ 2 e' 2 ^ dxdt 
Jo Jn Jo Jn 

<C\ Ff \ 2 (s6) 2 \q\ 2 e- 2sv dxdt+ Ff X\s9) 3 \p\ 2 e- 2sv dxdt j (4.23) 

\Jo JX(u,t,0) Jo Jx(ui,t,0) J 

for appropriate s > s 2 > s% and A > A2 > Ai. 
Step 2. Arrangements and conclusion. 

In this step we estimate the local integral of p appearing in (4.23) by a local integral of q and 
some small order terms. Finally, using semigroup theory, we finish the proof of Proposition 4.2. 
The main result of this step is the following. 

Lemma 4.6. There exist some numbers A2 > Ai, s 2 > Si and C 2 > such that for all A > A2, 
all s > s 2 and all (po,qo) £ L 2 {£1) 2 , the corresponding solution (p,q) of system (4.6)-(4.10) 
fulfills 

Ff [(s6)- 1 {\Ap\ 2 + \p t \ 2 ) + \ 2 (s9)\X7p\ 2 + \ 4 (s6f\p\ 2 ]e- 2sv dxdt+ Ft (A 2 s9)\q\ 2 "e~ 2s * 'dxdt 
Jo Jn Jo Jn 

< C 2 f I \ s (s9) 7 e- 2sip \q\ 2 dxdt. (4.24) 

Jo Jx{u),t,0) 

Proof of Lemma 4-6- In order to prove Lemma 4.6, we just need to estimate the p appearing in 
the right-hand side of (4.23 ). For that, we introduce the function 

C(M) := £(X(x,0,t)), (4.25) 
where £ is a cut-off function satisfying 

£eC£>), (4.26) 

< f (a?) < 1, x £ R N , (4.27) 

£(x) = 1, x 6 wi. (4.28) 

We have that 

Ff X\s9f\p\ 2 e- 2sip dxdt< Ff (\ A (s9f\p\ 2 e- 2sip dxdt (4.29) 
Jo Jx(uii,t,a) Jo Jn 

and we use (4.7) to write 



Ff (\ 4 (s9) 3 \p\ 2 e- 2sip dxdt = Ff (\\s9fpqe- 2sv dxdt 
Jo Jn Jo Jn 



+ f f (\ 4 (s9fp(-q t )e- 2sv dxdt =: M 1 + M 2 . (4.30) 
Jo Jn 
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It remains to estimate M\ and M^. Using Cauchy-Schwarz inequality and (4.26)-(4.27), we 
have, for every e > 0, 

|Mi|<e [ T [ \ 4 (s9f\p\ 2 e- 2sv dxdt + ^ [ T [ \ A (s6) 3 \q\ 2 e - 2 ^dxdt. (4.31) 



/o Jn 4e j Q jx(u,t,o) 

On the other hand, integrating by parts with respect to t in M 2 yields 

M 2 = f f C\ A {sefp t qe- 2stp dxdt+ [ [ (\ 4 (3s 3 6 2 6 t - 2s i ^> t e 3 )pqe' 2sLp dxdt 
Jo Jn Jo Jn 

- £j V£(X(x, 0, t)) ■ 0, t),t, 0)f(x, t)\\s6fpqe- 2s *dxdt 

= : m\ + m| — m|. 

For M2, we notice that for every e > 0, 

I M 2 X I < e f T f (s6)~ l \p t \ 2 e- 2sv dxdt + ±- A 8 (sfl) 7 |g|V 2s ^xcii. (4.32) 



Since + |(/?t| < CX8 , we infer that 

rT 



|M?I < 



C / f (s\\6) 5 \pq\e- 2s ^dxdt 
Jo Jn 

£s * Jo Jx(ui,t,0) 



< e f f \ 4 (s6f\p\ 2 e- 2sip dxdt + -^ f f \\s6) 7 \q\ 2 e- 2 ^dxdt. (4.33) 
Jo Jn £s Jo Jx(u,tjo) 



Finally, Mf is estimated as M±: 



4f m3u2 -2 S¥3 , , C f T l \i| C .av*u2,--'-v, 



\M%\<e \\s6y\p\ z e- zs *dxdt+- / A>0)%| 2 e -^cte<ft. (4.34) 

Gathering together (4.23) and (4.29)-(4.34) and taking e small enough, we obtain (4.24). 
Now we finish the proof of the observability inequality (4.11). 

Pick any (po,qo) G L 2 (Q) 2 , and denote by (p,q) the solution of (4.6)-(4.10). Note that 
p G C([0,T];L 2 (fl)) nL 2 (0,T;H%(n)) and that q £ H 1 ^, T; L 2 (f2)). Using classical semigroup 
estimates, one derives at once (4.11) from (4.24). 

□ 

5. Final comments 

• Another decomposition 

As commented in the introduction, there is another splitting of the operator L = 
df — A — Adt + dt, given by 

L = (d t -A)(d t + id). 

Thus, letting 

v(x,t) = y(x,t) +yt(x,t), 
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we see that (1.1) may be written as 

v t -Av = l u(t )h+(l-b(x))(v-y), (5.1) 
yt + y = v, (5.2) 

which is a coupled system of a parabolic equation (5.1) and an ODE (5.2). 

This splitting was used to prove Theorem 1.1 with less assumptions on the trajectories 
(see Section2). 

The control term h acts directly in the heat equation and indirectly in the ODE 
through the coupling term v. The problem can be treated directy as such, with requires 
further work at the level of the dual observability problem since both components of the 
adjoint system will be needed to be observed by partial measurements only on one of 
its components. The problem can also be addressed incorporating in (5.2) an additional 
auxiliary control acting directly in the ODE. This leads to the system 

v t -Av = l w(t) /i+(l-6(x))(«-y) (5.3) 

yt + y = l u (t)k + v, (5.4) 

where (v, y) G L 2 (tt) 2 is the state function to be controlled, and (h, k) G L 2 (0, T; L 2 (0) 2 ) 
is the control input. 

Once the controllability of this system is proved, when going back to the original 
viscoelasticity equation, one gets 

y tt -Ay- Ay t + b{x)y t = l w(t) [h-(l- b{x))k] + (8 t - A)[l u{t) k]. (5.5) 

But, then, the second control l w mfe enters under the action of the heat operator. It is 
then necessary to ensure that the control k is smooth enough and, furthermore, to replace 
in (5.4) the cut-off function luti) by a regularized version. These are technicalities that 
can be overcame with further work. To be more precise, the control in (5.5) takes the 
form l X (u,tfi){x)h with h G L 2 (0, T; L 2 (n)), provided that both h, k G L 2 (0, T; L 2 (£l)) 
and 

k G ^(0, T; L 2 (Q)) n L 2 (0, T; H 2 (Q)). 

Therefore special attention has to be paid to obtain smooth controls for the transport 
equation (see Section 2). 
• Manifolds without boundary 

The lack of propagation properties of the ODE (1.10) in the space variable requires 
the control to move in time. As we mentioned in the introduction, through a suitable 
change of variables, this is equivalent to keeping the support of the control fixed but 
replacing the ODE by a transport equation. Obviously, attention has to be paid to 
the Dirichlet boundary conditions when performing this change of variables. Of course, 
this is no longer an issue when the model is considered in a smooth manifold without 
boundary. As an example of such a situation we consider the periodic case in the torus 

xeT N := R N /Z N . (5.6) 

For a moving control with a constant velocity u(t) = {x — ct; x G w}, c G \ {0}, 
system (1.9)-(1.10) can be put in the form of a coupled system of parabolic-hyperbolic 
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equations 



vt — Av — c • Vv + (b(x + ct) — l)v 
wt — c ■ Vu> + w 



w 



l u (x)h + (b(x + ct) - l)v 



(5.7) 
(5.8) 



by letting 



v(x, t) 
w(x, t) 
h(x, t) 



y(x + ct,t), 
z(x + ct, t), 
h(x + ct, t). 



(5.9) 
(5.10) 
(5.11) 
(5.12) 



The system is now constituted by the coupling between a heat and a transport equation 
with control h with fixed support. Once more, the problem now can be treated by means 
of the classical duality principle between the controllability problem and the observability 
property of the adjoint system. The later was solved in [25] in 1 — d using Fourier analysis 
techniques and in this paper we do it using Carleman inequalities. 

Note that the Carleman approach developed in this paper cannot be applied as it 
is to the periodic case. Consider for instance the case of the torus T. A weight tp S 
C°°(T x (0, T)) as in Lemma 4.3 does not exist, because of the periodicity in x (see 
Figure 6.) However, it is well known that the periodic case can be deduced from both 
the Dirichlet case and the Neumann case (using classical extensions by reflection, see 
e.g. [ ]). Even if the Neumann case was not considered in this paper, it is likely that 
it could be treated in much the same way as we did for the Dirichlet case. 



The main result of this paper concerns the controllability of a coupled system consisting on 
a heat equation and an ODE. By addressing the dual problem of observability and making the 
controller /observer move in time, this ends being very close to the problem of observability of 
a coupled system of a heat equation and a first order transport equation. The techniques we 
have developed here are inspired in the work [1] where the key point was to use the same weight 
function for the Carleman inequality in both the heat and the transport equation. 

The system under consideration, coupling a heat and a hyperbolic equation, is close to that of 
thermoelasticity that was considered in [19]. But, there, the problem was only dealt with in the 
case of manifolds without boundary, by means of spectral decomposition techniques allowing to 
decouple the system into the parabolic and the hyperbolic components. As far as we know, a 
complete analysis of the system of thermoelasticity using Carleman type inequalities seems to 
be not developed so far. 

The structure of the parabolic-transport system we consider is also, in some sense, similar 
to the one considered in [9] for the 1 — d compressible Navier-Stokes equation although, in the 
latter, the system is of nonlinear nature requiring significant extra analysis beyond the linearized 
model. 

Our analysis is also related to recent works on the control of parabolic equations with memory 
terms as for instance in [13]. Note that the system (1.9)-(1.10) in the particular case 6=1 and 



6. Open problems and further questions 
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z(0) = 0, in the absence of the control h and the addition of a control of the form l^mk in the 
first equation, can be written as an integro-differential equation 

yt-Ay + (b-l)[y- I e^yix, s)ds] = l u{t) k (6.1) 
J o 

This system is closely related to the one considered in [13]. There it is shown that the system 
lacks to be null controllable. This is in agreement with our results that, in the particular case 
under consideration, show also that a moving control could bypass this limitation. It would 
be interesting to analyze to which extent this idea of controlling by moving the support of the 
control can be of use for more general parabolic equations with memory terms. 

In this paper we have shown the null controllability of a linear system which consists of a 
parabolic equation and an ordinary differential equation that arise from the identification of 
the parabolic and hyperbolic parts of system (1.1)-(1.2). Besides, coupled systems consisting of 
parabolic equations and ode's are important since they appear in biological models of chemotaxis 
or interactions between cellular process and diffusing growth factors (see [14], [24], [27] and 
references therein). Systems governing these phenomena are, in general, non linear and have 
the form 



u t = f(u,v), (6.2) 
v t = DAv + g(u,v), (6.3) 

where v and u are vectors, D is a diagonal matrix with positive coefficients and / and g are real 
functions. 

Other area where coupled parabolic-ode systems play a major role is electrocardiology (see [3], 
7], [33] and references therein). Here the cardiac activity is described by the bidomain model, 
which consists of a system of two degenerate parabolic reaction-diffusion equations, representing 
the intra and extracellular potential in the cardiac muscle, coupled with a system of ordinary 
differential equations representing the ionic currents flowing through the cellular membrane. 



The bidomain model is given by 

xC m v t - div (AVuj) + xhon{v, w) = r app , (6.4) 

- xCmVt - dlV (D e Vu e ) - xhon(v, V)) = -Iapp, (6.5) 

w t - R{v,w) = 0, (6.6) 



where u% and u e are the intra and extracellular potentials, v is the transmembrane potential, x 
is the ratio of membrane area per tissue volume, C m is the surface capacitance of the membrane, 
lion is the ionic current, I l dp P is an applied current and Di^ e are conductivity tensors. 

Concerning controllability of coupled parabolic-ode systems, just a few results for some par- 
ticular systems are known (see [8] and [23] for the controllability of a simplified one- dimensional 
model for the motion of a rigid body in a viscous fluid). We believe that ideas presented in 
this paper can be used for the study of the controllability for other systems of parabolic-ode 
equations, such as (6.2)-(6.3) and (6.4)-(6.6). 
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Appendix A. Proof of Lemma 4.3 

Proof. Pick any 5 < min(ii,T — t 2 ,T/2). We search ip (see Figure 6) in the form 

iP(x,t) = ijj 1 {x,t) + C 2 Mx,t) + C 3 (A.l) 

where, roughly, Vr fulfills (4.13), ip 2 fulfills (4.14)-(4.16) together with Vip 2 = outside X(co u t, 0), 
and C 2 ,Cs are (large enough) positive constants. 
Step 1. Construction of -0i- 

Let T G C°°([0,T\;R N ) be as in (1.13), and let e > be such that 

B(T(t),3e) C X(u ,t,0)nn, t € [0,T\. 
We choose a vector field / G C°°(R N x [0, T];M N ) such that 

f (T f\ = i r(*) ifte[0,T], xeB(T(t),e), 

J{ , ) \0 if t€ [0,T], sel w \B(r(t),2£). 

Let X denote the flow associated with /; that is, X solves 
dX 

— {x,t,t ) = f(X(x,t,t ),t), (x,t,t )eR N x [0,T] 2 , 

X{x, t , t ) = x, (x, t ) €R N x [0, T\. 

Note that 

A(y + r(0),t,0) = y + T(t) if (y, t) G £?(0,e) x [0, T], 

i"(x, t, t ) = ^ if dist (x,9Q) < e, (Mo) G [0,T] 2 . 

By a well-known result (see [16, Lemma 1.2]), there exists a function ip G C°°(f2) such that 

tjj(x) > if x GO,; 

i>(x) = if x G <9ft; 

Vip{x) ^ if x G O \ B(T(0), e). 

Actually, the function ^ given in [1G] is only of class C 2 , but the regularity C°° can be obtained 
by mollification with a partition of unity (see e.g. [29, Lemma 4.2]). Let us set 

Mx,t) = ip(X(x,0,t)). 

Then G C°°(Ti x [0, T\) and it fulfills 

^i(x,i)>0 if (x,t) G n x [0,T], (A.2) 

V>i(a?,t) = if (i,i)6 9Qx[0,T], (A.3) 
ay 

VA(aj,i) = VV>(X(x,0,t)) — (x,0,t) / if x G fi \ X(w , f, 0). (A.4) 

ox 

For (A.4), we notice that if we write x = X(x, t, 0), then x = A(x, 0, t) hence 

V^(A(x,0,t)) = V^(x) / 

if x S(r(0),e), which is equivalent to x -B(r(t),e). The last condition is satisfied when 
x X(u ,t,0). 
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Step 2. Construction of ip 2 . 

From (1.15), (1.16), and (1.17), we can pick two curves 71 G C°([0, i 2 ); fi) and 72 G C°((ii, T];Q) 
such that 



71 (t) 0X(w o ,*,O), 0<t<t 2 , 
72 (t) ^1M0), t x <t<T. 

We infer from (1.17) that for any f G (ii, f 2 ), 7i(0 and 72 (i) do not belong to the same connected 
component of O \ X(coo,t, 0). Let fii(i) (resp. 2 (i)) denote the connected component of 71 (i) 
(resp. 72 (t)) for < t < t 2 (resp. for t\ <t <T). Clearly 

( fii(t), if0<t<ti, 

n\X(wo,t,0) = < ni(i)Ufl 2 (i), ifti<*<t 2j 

i n 2 (t), if t 2 < t < r. 

Set n x (t) = for i G [t 2 ,T], and ft 2 (i) = for t G [0,ti]. Let ip 2 G C°°(n x [0,T]) with 
ih{x,t) = t (l Ul(t) (x) - l Q2(t) (x)) for t G [0, T], x E Q\ X(coi, t, 0), 

^ = for (x,t) G dfi X [0,T1. 
an 

Such a function ijj 2 exists, since by (4.12) 

inf dist (Q 1 {t)\X(u 1 ,t,0),n 2 (t)\X(uj 1 ,t,0)) > 0. 

Then 

chh _ ( 1 if 0<t<t 2 and x G fii(i) \ X(a; x , 0), 
dt ~ \ -1 if h<t<T and x G fi 2 (i) \ X(w 1; i, 0). 

and 

Vifj 2 (x,t) = if xGQ\X(wi,i,0). 

Note that (4.14)-(4.16) are satisfied for ip 2 . Note also that for any pair (ri,r 2 ) with < Ti < 
t 2 < T, the set 

K TUT2 := {(x,t) G R N+1 ; n < t < t 2 , x G Vl \ X(u u t, 0)} 

is compact. 

Step 3. Construction of if). 

Let ^ be defined as in (A.l), with C 2 > and C3 to be determined. Then (4.13) and (4.17) are 
satisfied. We pick C 2 large enough for (4.14)-(4.16) to be satisfied. Finally, (4.18) is satisfied for 
C3 large enough. □ 

Appendix B. Proof of Lemma 4.4 

Proof of Lemma 4-4- The method of the proof is widely inspired from [12], and the compu- 
tations are presented as in [29, Proof of Proposition 4.3]. 
Let v = e'^p and P = d t + A. Then 



e- a *Pp = e- s<p P(e Sip v) = P s v + P a v 
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where 

P s v = Av + (stpt + s 2 \Vip\ 2 )v, (B.l) 

P a v = Vt + 2sV<p ■ Vv + s(Aif)v (B.2) 

denote the (formal) self- adjoint and skew-adjoint parts of e~ S(p P(e sip -), respectively. It follows 
that 

\\e~ sip Pp\ | 2 = | \P s v\ | 2 + ||P a t>|| 2 + 2{P s v, P a v) (B.3) 

where (f,g) = f^fgdxdt, \\f\\ 2 = (/,/). In the sequel, J Q T J n f(x,t)dxdt is denoted fff for 
the sake of shortness. We have 

(P s v, P a v) = (Av, v t ) + (Av, 2sVif ■ Vv) + (Av, s(A(p)v) + (sip t v + s 2 |V^| 2 -y, v t ) 

+ {stp t v + s 2 \V(p\ 2 v, 2sV(p ■ Vv) + (sip t v + s 2 \Vip\ 2 v , s(A<p)v) =: h + h + h + h + h + h- 

(B.4) 

First, observe that 

h = - jj Vv-Vv t = 0. (B.5) 
Using the convention of repeated indices and denoting di = d/dxi, we obtain that 
h = 2s JJd 2 vdi(pdiV 

= —2s / / djv(djdifdiV + dupdjdiv) + 2s / (djv)rijdi(pdivda. 
J J Jo Jdn 

Since v = for (x,t) £ d£l x (0,T), Vv = (dv/dn)n, so that Vip ■ Vv = (d<p / dn)(dv / dn) and 

n(djv)rijdi(pdiV da = / / (dip/dn)\dv/dn\ 2 da. 
m Jo Jon 

It follows that 

I2 = —2s / / djdiipdjvdiV — s / / di(pdi(\djv\ 2 ) + 2s / / (d(p/dn)\dv/dn\ 2 da 
J J J J Jo Jon 

= -2s [ [ djdifdjvdiv + s [ [ A^\Vv\ 2 + s [ [ (d^/dn)\dv/dn\ 2 da (B.6) 
J J J J Jo Jen 

On the other hand, integrations by parts in x yields 

I 3 = -s JJvv- (vV(A<p) + (A^)Vv) = s jj A 2 v?^- - s jj A^|Vv| 2 (B.7) 
and integration by parts with respect to t gives 

h = - ff{s^ t + s 2 d t \V^\ 2 )^- 
Integrating by parts with respect to x in ^5 yields 

h = ~ jj s 2 V ■ (vtVtfM 2 - jj 's 3 V- (iVvfV^M 2 . (B.8) 
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Gathering (B.4)-(B.8), we infer that 
2{P s v,P a v) = -4s If d j d i (pd j vd i v + 2s I I (dip / dn)\dv / dn\ 2 da 

J J JO Jdtt 

+ J J |w| 2 [s(AV - <Pu) - 2s 2 ft|V^| 2 - 2s 3 V^ ■ V|V(^| 2 ]. 
Consequently, (B.3) may be rewritten 

\\e- Sip Pp\\ 2 = \\P s v\\ 2 + \\P a v\\ 2 -As II djdupdjvdiv + 2s I I (dip/dn)\dv/dn\ 2 da 

J J JO J dfi 

+ J J | V | 2 [ S (A 2 ^-^)-2 S 2 a 4 |V^| 2 -2s 3 V^-V|V^| 2 ]. 

Claim 1. There exist some numbers Ai > 0, s\ > and A G (0, 1) such that for all A > Ai 
and all s > s±, 

\v\ 2 [s{A 2 ip - tp tt ) - 2s 2 d t \Vip\ 2 - 2s 3 V<p • V|V<^| 2 ] 

1\„3 / / / \M3i„,|2 v. A\ n 3 I I f\a\3\„,\2 



+ A- L \s 6 / (X0y\v\ z > AXs 3 / / (\6y\v\. (B.9) 

Jo Jx(uj 1 ,t,0) J J 

Proof of Claim 1. Easy computations show that 

dicp = -Xg(t)e^d^, Ojdtif = -g(t)e^(X 2 d l ipd^ + Xdjd^) (B.10) 

and 

-V|Vv9| 2 • V99 = -2{djdnp)dnpdw = 2(A 5 e AV ') 3 (A|VV'| 4 + djd^d^d^). 
It follows from (4.13) that for A large enough, say A > Ai, we have that 

-V\V<p\ 2 -Vcp > AX(X6) 3 , t G [0,7], x G U\X(ui,t,0) (B.ll) 

|V|Vc^| 2 - Vv?| < A^XiXOf, t G [0, 7], x G X(ui\,t, 0) (B.12) 

for some constant A G (0, 1). According to (4.18), we have for some constant C > 

\A 2 <p\ + \ip tt \ + |d t |V^| 2 | < CX(X9) 3 , t G [0,T], x G H. 

Therefore, we infer that for s large enough, say s > si, and for all A > Ai we have that 

s(A 2 cp-ip tt ) -2s 2 d t \Vip\ 2 -2s 3 Vip- V\Vip\ 2 > AXs 3 (X9) 3 , t G [0,7], x G n\X(wi,t,0) 

|s(AV -Vtt) -2s 2 d t |Vv?| 2 -2s 3 V<^- V|V<^| 2 | < 3A- 1 Xs 3 {X6) 3 , t G [0,T], x G X(ui,t,0). 

This gives (B.9) with a possibly decreased value of A. □ 
Thus, using the fact that dip/dn > on dQ by (4.17), we conclude that 

||i>|| 2 + \\P a v\\ 2 + AXs 3 JJ (A#) 3 M 2 

<\\e- slp Pp\\ 2 + 4s II djdiipdjvdiV + A^Xs 3 I I (A#) 3 M 2 . (B.13) 

J J Jo JX(ui,tJ0) 
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Claim 2. There exist some numbers A2 > Ai and S2 > s% such that for all A > A2 and all s > S2, 
As JJ (X9)\Vv\ 2 + Xs' 1 JJ (X9)- 1 \Av\ 2 < C (s- l \\P s v\\ 2 + \s z JJ (A6>) 3 |^| 2 ^ . (B.14) 
Proof of Claim 2. By (B.l), we have 

s- 1 JJ (X9)~ 1 \Av\ 2 = s- 1 Jjixey^Psv-siptv-s^vcpfv] 2 

< cs- 1 J f {\e)-\\p 8 v\ 2 + s 2 M>| 2 + s\\6) 4 \v\ 2 ) 

< C (p^f*- + s 3 j "j \\ef\v\ 2 ) (B-15) 

provided that s and A are large enough, where we used (4.18) in the last line to bound ipt- On 
the other hand, 

As JJ {X9)\Vv\ 2 = Xs{JJ (X9)(-Av)v - JJ (V(A0) • Vv)v] 

< A II (xer^v\ 2 + (\ef\v\ 2 + ±[[ A(xe)\ v \ 2 



2s JJ 2 JJ x J ' ' 2 _ 

< c (s-^PsvW 2 + Xs 3 JJ(xe) 3 \v\ 2 j (B.16) 

by (B.15), provided that s > S2 > si and A > A2 > Ai. Then (B.14) follows from (B.15)- 
(B.16). □ 
We infer from (B.13)-(B.14) that 

\\P a v\\ 2 + Xs IJ(Xe)\Vv\ 2 + Xs- 1 JJ(X9)- 1 \Av\ 2 + Xs 3 JJ(X9) 3 \v\ 2 

<C[\\e- sip Pp\\ 2 + 4:S II d j dnpd j vd i v + A- l Xs 3 I I (X6) 3 \v\ 2 ) . (B.17) 
y J J Jo Jx(ui,t,o) J 

By (B.10), 

s JJ djdipdjvdiv <-sX JJ g(t)e H 'djd^djvdiv <Cs JJ (X9)\Vi 
Therefore, for A large enough and s > S2, 

\\P a v\\ 2 + Xs 3 II (X6) 3 \v\ 2 + Xs II (X6)\Vv\ 2 + As" 1 / / (X9)- 1 \Av l2 



>v\ 2 . 



<C\\\e- SLp Pp\\ 2 + Xs 3 T I {X6) 3 \v\ 2 ] . (B.18) 

V Jo ix(wi,i,0) / 
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Using (B.2) and (B.18), we see that for A large enough and s > S2 

A*- 1 !!(\ey l \v t \ 2 < cxs- 1 !! (\er l {\p a v\ 2 + s 2 \v^\ 2 \vv\ 2 + s 2 \^\ 2 \ v \ 2 ) 



T 

< C I \\e-^Pp\\ 2 + Xs 3 I / (A0) 3 M 2 
y Jo Jx(ui,t,o) J 

Hence, there exists some number A3 > A2 such that for all A > A3 and all s > S2, we have 
As 3 f [ (\6f\v\ 2 + As ft (X0)\Vv\ 2 + As" 1 f [ (Xe^QAvf + \v t \ 2 ) 



<C[\\e- sip Pp\\ 2 + Xs 3 [ T [ (X9) 3 \v\ 2 ) . (B.19) 
y Jo Jx(ui,t,o) J 



Replacing v by e sip p in (B.19) gives at once (4.19). The proof of Lemma 4.4 is complete. □ 

Appendix C. Proof of Lemma 4.5 

Proof of Lemma 4-5. The proof is divided into three parts corresponding to the estimates for 
t € [0, 5} , for t £ [5, T - 5} and for t G [T-5,T\. The estimates for t 6 [0, 5] and for i e [T - S, T] 
being similar, we shall prove only the first ones. 

Let v = e'^q. Then 

e-Vqt = e-^(e s ^v) t = s<p t v + v t =: P s v + P a v. (C.l) 

Claim 3. 

" X(s6) 2 \v\ 2 dxdt < C ( / / X-^e-^qtfdxdt 



oJn \JoJn 



+ [ [(l-() 2 (s6)\v\ 2 ] lt=s dx+ f S f X{s6) 2 \v\ 2 dxdt j , (C.2) 

if] JoJx(uj,t,0) J 



where £ is the function introduced in (4.25). 

To prove the claim, we compute in several ways 

I:= [ [ {e-^qt^l-QfsOvdxdt. 
Jo Jo. 

We split / into 

1= [ [ {P s v){l-C) 2 s0vdxdt+ ( ( {P a v)(l-C) 2 s9vdxdt=:h + I 2 . 
JoJn Jo Jo. 

Then 

(l-() 2 s 2 (p t 6v 2 dxdt 

[ 5 '( e |A||^IUoc _ e A^ _ g x^ te ^](i - () 2 s 2 ge x ^v 2 dxdt. 

oJn 
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On the other hand 



h = f /(I -C) 2 (sge Xi, vv t )dxdt 



' f [(l-0 2 s6\v\\ =5 dx- I* I s[g'e^ + g\^]{l-C) 2V *-dxdt 



' (1 - C)V£(X(x, 0, 0) • (^)'\x( X , 0, t), *, 0)/O, t)s9v 2 dxdt, 



dx 



where we used the fact that [#|-u| 2 ]| <=0 = 0. Clearly, since 9 > 1, for s > 1 

<5 /• ay 

'1 - C)V£(X(x, 0, t)) ■ (—)~ (X(x, 0, i), t, 0)f(x, t)s9v 2 dxdt\ 

ox ' 

<C ( ( (s9) 2 \v\ 2 dxdt. 

On the other hand, using (4.15), we see that there exist some constants C > and s\ > sq such 
that for all s > s\ and all A > Ao > 0, it holds 

gXip t e Xi, (s 2 ge x ^ + |) > CA(s6>) 2 , t G (0, 5), x G H \ X(u} 1} t, 0) 

- 5 '(i) ((el A|W ^° - e x ^)s 2 ge^ - V^) > t G (0, <J), x G H \ X(wi, t, 0). 

It follows that for some positive constant C > C 



1 



5 



C7 / / A(s#)>| 2 dxcft < -I+- / \{1-Cf{s9)\v\\ =s dx + C' / / X(s9) 2 \v\ 2 dxdt. (C.3) 

Join 2 Jj) JoJx{u,t,o) 

Finally, by Cauchy-Schwarz inequality, we have for any k > 

|/| < (4k)- 1 I ! \e~ sv q t \ 2 dxdt + k I I (s9) 2 \v\ 2 dxdt. (C.4) 
Jo Jn Jo Jn 

Combining (C.3) with (C.4) gives (C.2) for n/X > small enough. Therefore, Claim 3 is 
proved. □ 
We can prove in the same way the following estimate for t £ [T — 5,T]: 

T 



T—SJfl 



X(s9) 2 \v\ 2 dxdt <c( I [ X- l \e- sv q t \ 2 dxdt 
\Jt~5 Jn 

+ f [(l-() 2 (s9)\v\\ =T _ 5 dx+ [ T [ X{s9) 2 \v\ 2 dxdt 
Jn JtsJ x(u;,t,o) 



(C.5) 
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Let us now consider the estimate for t £ [5, T — 5] . 
Claim 4. 

' \ 2 (s6)\v\ 2 dxdt + I [(l-0\Xse)\v\\ =5 dx+ I [(l-0 2 (\s6)\v\\ =T „ s dx 
n Jn Jn 

f [ \e- Sip qt\ 2 dxdt+ [ [ \ 2 (s6)\v\ 2 dxdt) . (C.6) 

Js Jn Js Jx(u>,t,o) J 

|| • || and (., .) denoting here the Euclidean norm and scalar product in L 2 ($7 x (5,T — 5)), we 
have that 

He-^H 2 > ||(1 - C)(W + vt)\\ 2 > 2((1 - C)W, (1 " CM)- (C7) 
Next, we compute 



((i - (i - CM) = y (i - CYs Vt -dx 

T-S r dX 



TS ,-T-S 
t=S ~2 



(1 - C) 2 ^|w| 2 rfxdt 



(l-C)Ve(X(x,0,t))-(— ) 1 (X(x,0,t),t,0)/(x } t)s^i; 2 da:dt. (C.8) 

Since g(t) = 1 for i € [6,T — 5], we have that cpt = — A^je A ^. From (4.15)-(4.16), we infer that 
sip t (x,T-5) > C\se x ^ X(oj\, T — 5,0), 

-s(f t (x,6) > CXse x ^ x \ X(uj\, 5, 0). 
Therefore, using (4.28), 



2 V 

(1 - C) stpt—dx 



T-S 



>c( I [(l-C) 2 (Xs0)\v\\ =s dx+ [ [(l-() 2 (\s8)\v\ 2 ] lt=T _ s dx 
VJn Jn 



(C.9) 

Next, with tp u = -{(\ip t ) 2 + Xip tt }e x ^ and (4.14), we obtain for A > Ai > A 

- S - f [ (1 - Q 2 iptt\v \ 2 dxdt > C [ [ (1 - C) 2 A 2 ^|u| 2 dxdt. (CIO) 

2 J«5 Jf7 J<5 Jf7 

Finally 

y (1 - C)V£(X(x, 0, t)) • (— 0, i), t, 0)f(x, t)s^v 2 dxdt 

<C I [ \s6\v\ 2 . (C.ll) 

Claim 4 follows from (C.7)-(C.ll). 

We infer from (C.2), (C.5) and (C.6) that for some constants Ai > Ao, s\ > sq and C\ > we 
have for all A > Ai and all s > s± 

/ ( \ 2 {se)\v\ 2 dxdt<Ci[ I' [ \e~ slp qt\ 2 dxdt+ [ [ \ 2 {s9) 2 \v\ 2 dxdt J . (C.12) 
io Jn yJo Jn Jo Jx(u,t,o) J 

Replacing u by e~ sip q in (C.12) gives at once (4.20). The proof of Lemma 4.5 is complete. □ 
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